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Abstract 

Possible effects of the color screened confinement potential are investigated. Color 
screened linear potential with a large string tension T = (0.26 — 0.32)GeV 2 is suggested by 
a study of the cc and bb spectra. The -0(4160) and 0(4415) are respectively assigned as the 
^(45)-dominated and the ip(5S) cc states. Satisfactory results for the masses and leptonic 
widths (with QCD radiative corrections) of cc and bb states are obtained. 
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The string tension is a most fundamental physical quantity in quark confinement. 
It is argued based on a rotating string picture that the string tension T is related to 
the Regge slope a' by M 



T= 2W = °- 18GeV2 > W 
with the experimental value for a = 0.9(GeV)~ 2 . This result is supported phe- 
nomenologically by the heavy quarkonium spectrum when identifying the string ten- 
sion with the slope of the linear confinement potential between a heavy quark- 
antiquark (QQ) pair! 2 !. Lattice QCD calculations for the string tension is not con- 
clusive, because one needs to estimate the lattice scale A L in physical units, even if 
one has obtained from lattice calculations the value for c which is related to T via 

In ref. [3] the string tension is estimated (though with some theoretical uncertainty) 
to be 

T = (0.33t° f 3 )GeV 2 (3) 

with typical values 

c = (7.5 ± 0.5) x 10- 2 , = (200^°)iWW, (4) 

although this estimate has some uncertainties, e.g., the value of c is calculated in 
the quenched approximation while the experimental value of is obtained with 
dynamical fermions (the more recent data give AgL = (195 + 65 - 50)MeV W). The 
central value in Eq.(3) is significantly larger than that given in Eq.(l), but they may 
be still consistent with each other when the large error involved is reduced. 

On the other hand, while the potential model, which assumes a linear confine- 
ment potential plus a one-gluon exchange potential, is generally successful for cc and 



bb spectroscopy, some problems may remain. One of them is concerned with the 
assignment of ■0(4160) and ■0(4415). With the linear confinement potential, -0(4160) 
and -0(4415) are usually assigned as the 2D and 4S* states respectively. However, 
experimentally the 0(4160) has a quite large leptonic width M T ee = 0.77 ± 0.23keV, 
comparable to that for the 3S state 0(4040), but in the nonrelativistic limit the 2D 
state will be forbidden to decay to e + e~. Neither the S — D mixing nor the coupled 
channel models can consistently solve this problem' 5 '. Moreover, for 0(4415), with 
the linear potential the leptonic width of AS state is usually predicted to be larger by 
more than a factor of two than the observed value of -0(4415), hence the assignment 
of -0(4415) as AS state is also problematic. The way to solve this problem is probably 
to assign the 0(4160) and -0(4415) as the AS and 5S states respectively, but this must 
require the linear potential to be softened at large distances. In addition, with the 
linearly rising potential the calculated masses and leptonic widths for highly excited 
bb states (e.g. the 6S state) usually are also larger than their observed values. It is 
clear that for a pure linear potential all the wave functions (S-wave) at origion take 
the same value. Therefore, if the QQ potential keeps linearly rising at large distances 
the leptonic widths for highly excited states will gradually approach to a constant 
value when the linear potential becomes dominant over the short ranged Coulomb 
potential. However, the experimental data for both cc and bb do not show this ten- 
dency at all. On the contrary, as argued above the data stronly indicate a softening 
effect on the linear potential. 

Indeed, it is expected (e.g. in the color flux tube picture for confinement) that at 
large distances the creation of a light quark pair will screen the static color sources 
of QQ, and will therefore flatten the linear potential' 6 '. Some recent lattice QCD 
calculations with dynamical fermions seem to indicate that the color screening effects 
on the linear potential do exist at large distances' 7 ] . 

By taking into account the color screening effect, the QQ potential may be mod- 
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ified and take the form 

An 1 — p~^ r 

where the first term on the right hand side is the usual one gluon exchange Coulomb 
potential, and the second term is the screened confining potential with a screening 
parameter \i. The potential will keep linearly rising up to a distance r < [T x ( 
r < (1 — 2)fm for /i = (0.2 — 0.1)GeV ) and then gradually become flattened and 
eventually reach a constant value ^ . 

Since the presence of the Cornell potential l 8 ' , many improved potentials have been 
suggested. In particular, the running coupling constant a s (see e.g. ref. [9]) and one 
loop QCD radiative corrections ^ have been successfully incorporated into the QQ 
potential. On the other hand, however, the understanding of the confining potential 
is still poor, though some phenomenological confining potentials have been considered 
to improve the fit to the heavy quarkonium spectra including that of the higher lying 
states' 11 ' 12 !. Since color confinement is the most important part of dynamics in hadron 
physics, it is necessary to have further studies regarding the confining potential to- 
gether with the QQ spectroscopy. In this connection the screened confining potential 
expressed in (5), i.e., V sc (r) = Tr[ }~ e ^ ' ) may differ from many other phenomeno- 
logical potentials (e.g. various power law potentials, see refs. [11,12] and references 
therein) in the following respects: 

(i) It is indicated by some lattice QCD calculations with dynamical fermions (see, 
e.g., ref. [7]). 

(ii) It keeps linearly rising up to about one fm and then gradualy becomes a 
constant at large distances, and therefore it incorporates the large distance asymptotic 
behavior of color screening into the linear confinement in a natural manner. This is 
well motivated theoretically. 

(iii) As noted previously' 13 ^ , the inclusion of color screening is connected to the 



removal of the infrared divergences of the QQ interaction kernel in the momentum 
space. In fact, in momentum space the screened confining potential reads' 13 ^ 

Vsc(P) = --5 3 (P) + • (6) 

» 71 (P +/i 2 ) 2 

In view of the regularization of the linear potential in momentum space, the form of 
(6) and hence (5) seem to be quite natural and unique. The non-vanishing value of 
the cut off \i is expected to be related to the polarization of dynamical light quark 
pairs. 

Although the exact form of confinement interaction has not been analytically 
derived from the first principles of QCD, we believe that the screened confining po- 
tential V sc expressed in (5) should be a better candidate for describing confinement 
than many other potentials. Potential (5) may have phenomenological implications 
and it has been used in the study of heavy flavor mesons. ' 13 1. It will be interesting to 
have further phenomenological investications regarding this potential. 

In the following we will use potential (5) to calculate the cc and bb mass spectra, 
and then find the possible phenomenological values for the string tension T. As the 
first trial in a previous paper' 14 ! ; we used 

T = 0.21GeV 2 , a s = 0.51, fj, = O.llGeV, m c = lAGeV, (7) 

as inputted parameters to solve the nonrelativistic Schrodinger equation with po- 
tential (5). The obtained mass spectrum is satisfactory with -0(4160) and -0(4415) 
assigned as 0(45") and ip(5S) respectively. However, there are two problems for (7). 
The first one is that with the same value for T and /i and a smaller value for a s , we 
cannot find good result for the bb mass spectrum. The second is that the value of a s 
in (7) seems too large, not compatible with the present value of QCD scale parameter 
(see below). 
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We now find that with 



T = 0.32GeV 2 , a s = 0.306, fi = 0.156GeV, m c = 1.6GeV (8) 

for cc, and 

T = 0.32GeV 2 , a s = 0.275, fi = 0.l32GeV, m b = A.8GeV (9) 

for bb, as the the parameters in potential (5) to solve the Schrodinger equation, good 
results for both cc and bb can be obtained. The calculated masses and leptonic widths 
for cc states are shown in Table 1, and for bb states in Table 2. The experimental data 
are given by the Particle Data Group' 4 '. The leptonic widths are calculated using the 
nonrelativistic expressions without QCD radiative corrections ( TgJ and with QCD 
radiative corrections ( r ee )(see, e.g., refs. [8,15,16]) 

Tee = 16 ™ 2 4^> (10) 

16 

r ee = r° e (i - —a s (m Q )l (11) 

where a s (mg) stands for the coupling constant at the QQ mass scale, and it can 
be determined in the time-like processes of heavy quarkonium decays. Here we use 

[16] 

a s (m c ) = 0.28 for cc and a s {m b ) = 0.19 for bb . These values of the running coupling 
constant are consistent with the QCD scale parameter A-^ 200MeV^ 16 l 

From Table 1 we can see that the ■0(4160) and 0(4415) are assigned as the AS 
and 5S states. The predicted leptonic widths for these two states are in excellent 
agreement with data, whereas in the usual potential models without color screening 
effects the ■0(4160) and -0(4415) are assigned as 2D and AS states and then -0(4160) 
would have zero leptonic width ( in the nonrelativistic limit ) and -0(4415) would 
have a leptonic width of, say l.lkeV^\ too large by more than a factor of two than 
the observed value (0.47 ± 0.10)/ceV. As for the mass of the AS state, the predicted 
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value is higher than ■0(4160) by lOOMeV, and this could be due to the neglect of 
S — D mixing and coupled channel effects. In any case, if 0(4415) is the 5S state, 
the 0(4160) must be a 4S'-dominated state with possibly some mixed components of 
2D and virtual charmed meson pairs. It might be interesting to note that in these 
assignments the cc would have an anomalous mass relation that m(AS) — m(3S) is 
smaller than m(5S) — m(AS). Exactly the same anomalous mass relation is also 
observed for the bb states^. These anomalous mass relations may imply that in the 
energy region just above thresholds of many opened channels (e.g., in 3.8— A.3GeV for 
cc) the masses of resonances can be significantly distorted. Of course, in explaining 
these difficulties there could be other possibilities such as the ccqq states [17] or ccg 
states' 18 ' . However, these states in general do not seem to have large enough leptonic 
widths to be the -0(4160), because their couplings to the photon are expected to be 
suppressed. In Table 1 the predicted mass for tp(3S) is now 4.03GeV, much closer to 
its observed value than A.llGeV^ predicted by usual potential models. Moreover, in 
Table 1 the predicted leptonic widths for ip(lS) and 0(25) also agree with data. 

From Table 2 we see that in general the calculated masses and leptonic widths for 
the bb states are also in good agreement with data. As a result, using potential (5) 
with parameters (8) and (9), the obtained cc and bb spectra are remarkably improved. 

We have also tried to fit the spin-averaged mass spectrum using potential (5). 
Here the spin-averaged masses for the S wave states mean the masses before hyperfine 
splittings, e.g., for cc m(lS) = \[3m{ J ftp) + m(rf c )]. For other S wave states, because 
of the lack of observed values for the 0~ mesons, we use calculated hyperfine splittings 
(see (12)) and observed 1~ meson masses to determine the spin-averaged masses. We 
find that with slightly adjusted parameters (e.g., a slightly smaller string tension and 
a slightly larger a s ) we can get good fit for the spin-averaged cc mass spectrum and 
leptonic widths. Again, the assignments of -0(4160) and -0(4415) as the AS and 5S cc 
states seem to require a screened confining potential with a large string tension. 
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We have also used a modified Coulomb potential (with a running coupling constant 
(x s (r)) and the screened confining potential to fit the heavy quarkonium spectra, and 
the obtained results are similar to that obtained with the fixed coupling constant a s . 
Namely, a large string tension with color screening is still needed if a good fit for 
the higher excited states is required. In another words, taking a running a s does not 
change the basic feature of our observation on the screened string tension, though a 
slightly smaller value, say T = (0.28 - 0.30)GeV 2 is found. 

These studies might indicate that the color screened quasi-confinement potential 
with a large string tension, say, T = (0.26 — 0.32)GeV 2 should be an interesting 
possibility. 

The following observations might be in order. 

(1) . In order to get better results for higher excited QQ states (e.g. ■0(4160) 
and -0(4415)), a screened confining potential plus a Coulomb potential (with fixed or 
running a s ) seem to work well. Whereas the unscreened linear potential give too large 
level spacings and leptonic widths. While a large a s with a normal string tension (e.g. 
as shown in (7)) is possible, a smaller a s , which is more consistent with the value 
of QCD scale parameter, with a larger string tension (e.g. as shown in (8) and (9)) 
seem to work better for both cc and bb states. The screening parameter \x is found to 
be (0.14 ± 0.03)GeV. This value is consistent with lattice QCD calculations' 7 '. 

(2) . In our calculation we have simply focused on the spin-independent solu- 
tions of the Schrodinger eqution including the mass spectra and the leptonic widths, 
and ignored the coupled channel effects and relativistic corrections. In fact, the 
coupled channel effects (see e.g. refs.[5,8]) and the relativistic corrections (see e.g. 
refs. [15,19,20,21,22]) within the linear potential model seem to be unable to solve the 
puzzle regarding -0(4160) and -0(4415), as well as some other highly excited states. 
For instance, in the linear confinement model of ref.[21], with relativistic corrections 
the masses are found to be (in units of MeV) 3097, 3527, 3681, 3846, 4108, and 4446, 
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for IS, IP, 2S, ID, 3S, and 4S* cc states respectively. We see that although relativis- 
tic corrections for cc are important ( the energy shifts due to relativistic corrections 
ranging from —61MeV to —219MeV from IS to 45* states), the energy spacings with 
relativistic corrections can be very similar to that obtained without relativistic cor- 
rections e.g. in the Cornell model' 8 ' (the model in ref. [21] will of course have different 
parameters from that in ref. [8]). This may imply that as far as the energy spacings 
are concerned the relativistic effects may be largely absorbed by the readjustment of 
potential parameters ( e.g., the value of string tension takes T = 0.22GeV 2 in ref. [21] 
while T = O.lSGel^ 2 in ref. [8] ) and therefore the relativistic effects appear to be 
small in practice. Although this result is seen specifically in the model of ref. [21], the 
conclusion here can be quite general and similar observations have also been made by 
other authors (see, e.g., ref. [22]). Hence the relativistic corrections with unscreened 
confining potential are expected to be not very helpful in solving the difficulties as- 
sociated with e.g. -^(4160) and -?/>(4415). So it is very likely that in order to improve 
the fit to the higher excited states the screened confining potential is still needed even 
with these coupled channel and relativistic effects taken into consideration. 

(3). We have tried to calculate the spin-dependent splittings of these heavy 
quarkonium states in a very simple version. If the spin-spin force is entirely due to 
the lowest order perturbative one-gluon exchange, the 0~ — l - meson mass splitting 
A will be given by 



A = ^?l*(0)| 2 . (12) 



327T« S 

Then for the J/tp and i] c , with a s = 0.306, m c = 1.6GeV as given in (8), and 
the Schrodinger wave functions obtained by using (8), we get a mass spliting A = 
llOMel^, slightly smaller than its experimental value (118 ± 2)MeV. As for the fine 
splittings the situation is more complicated, since the long-range nonperturbative 
forces may contribute. It is argued' 23 ' based on a consistent condition due to Lorentz 
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invariance that the confining potential should transform as a Lorentz scalar and there- 
fore induce a spin-orbit term which then compensates the short-ranged spin-orbit force 
caused by one-gluon exchange. On the other hand, in many studies regarding chiral 
symmetry breaking, in order to preserve chiral invariance the vector ( or at least 
the time component of a 4- vector) confining force has to be chosen ' 2 4. In the phe- 
nomenological studies of heavy quarkonium spectra (see, e.g., refs. [15,19,20,21,25,26]), 
though the scalar confining is favored a vector-scalar mixture for the confining poten- 
tial may work even better. For instance, the vector-scalar mixed confining potential 
in some models can give an excellent fit to the low lying cc and bb spectra' 25 '. It is 
also argued that the observed tiny mass difference ( about — 0.9MeV )W between the 
center of gravity of triplet IP and the singlet IP charmonium states may not nec- 
essarily mean the short-ranged perturbative hyperfine splitting is dominant because 
nonperturbative and other effects could be also important [25, 26]. Although there are 
uncertainties for the spin-dependent splittings and the Lorentz transformed structure 
of the confining potential especially the color screened quasi-confining potential, we 
believe it should be dominated by the scalar with possibly a small mixture of the 
vector component. We find for the P-wave mass splittings, if the screened confining 
potential is a pure scalar, then with (5), (8), and (9) the obtained splittings are too 
small. If it is a vector-scalar mixture with a weight factor being about 3:7 then a fairly 
good fit can be obtained. However, we would like to emphasize that to calculate the 
spin-dependent splittings the naive calculation given here with a Coulomb potential 
of fixed a s in (5) with (8) and (9) should not be a good one, and a more refined 
calculation with higher order perturbative corrections and nonperturbative effects is 
apparently better. Therefore our calculation for the spin-dependent splittings is not 
conclusive, and we will leave this to a more refined work. 

To conclude, by studying the heavy quarkonium spectra especially for the higher 
excited states, e.g. ■0(4160), ^(4415), and T(11020), we find some evidence for the 
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color screened confining potential, which is expected theoretically when the creation 
of dynamical light quark pairs at large distances is taken into consideration. A large 
string tension, say T = (0.26 — 0.32)GeV 2 is favored by an overall fit to the mass 
spectra and leptonic widths. The existing calculations with the coupled channel effects 
and relativistic corrections based on the unscreened linear potential model seem to 
unable to solve the difficulties associated with those higher excited states. Therefore 
the color screened linear confining potential with a large string tension should be an 
interesting possibility and deserve further investigations. 

This work was supported in part by the National Natural Science Foundation of 
China and the State Education Commission of China. 
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Table 1: Calculated masses and leptonic widths for charmonium states with the 
screened potential (5) and parameters (8), where T ee = T^ e (l — ^a s (m c )) with 
a s (m c ) = 0.28t 16 l 



States 


Mass(MeV) 


r° e (keV) 


r ee (keV) 


r-f-(keV) 


Candidate 


IS 


3097 


10.18 


5.34 


5.26 ±0.37 


^(3097) 


2S 


3686 


4.13 


2.17 


2.14 ±0.21 


^(3686) 


3S 


4033 


2.35 


1.23 


0.75 ±0.15 


^(4040) 


4S 


4262 


1.46 


0.77 


0.77 ±0.23 


^(4160) 


5S 


4415 


0.91 


0.48 


0.47 ±0.10 


^(4415) 


IP 


3526 








X(3526) c . . 9 


ID 


3805 








^(3770) 


2D 


4105 











Table 2: Calculated masses and leptonic widths for bottomonium states with 
the screened potential (5) and parameters (9), where T ee = r° e (l — |fa! s (m&)) w ith 
a s (m b ) = 0.191 16 !. 
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States 


Mass(MeV) 


r° e (keV) 


r ee (keV) 


r^(keV) 


Candidate 


IS 


9460 


1.94 


1.31 


1.32 ±0.03 


T (9460) 


2S 


10023 


0.90 


0.61 


0.58 ±0.10 


T (10023) 


3S 


10368 


0.62 


0.42 


0.47 ±0.06 


T (10355) 


4S 


10627 


0.47 


0.32 


0.24 ±0.05 


T (10580) 


5S 


10833 


0.37 


0.25 


0.31 ±0.07 


T (10860) 


6S 


11002 


0.30 


0.20 


0.13 ±0.03 


T (11020) 


IP 


9894 








X b (9900) c . o . 9 . 


2P 


10267 








Xb (10261) c . o . g . 


ID 


10152 










2D 


10451 
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